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abstract 


The  evolution  of  the  power  spectrum  of  surface 
gravity  waves  is  described  by  means  of  a  transport 
equation.  The  effects  of  a  slowly  varying,  prescribed 
ocean  current  and  nonlinear  wave-wave  interactions  are 
included.  A  definition  due  to  Wigner  of  a  localized 
power  spectrum  is  used  to  derive  the  transport  equation 
from  the  dynamical  equations  describing  surface  wave 


motion . 
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I.  Introduction 


In  this  paper  we  present  a  derivation  of  a  transoo:.  t 
equation  to  describe  the  evolution  of  the  power  spectrum  of 
surface  waves  in  the  deep  ocean.  Our  transport  equation  is 
obtained  directly  from  the  dynamical  description  of  surface 
gravity  waves  and  is  to  be  contrasted  with  more  empirical  formu¬ 
lations  (Harnett,  3968;  Thomson  and  West,  1973).  Included 
will  be  the  effect  of  a  slowly  varying,  prescribed  current 
and  the  influence  of  nonlinear  wave-wave  interactions.  The 
dynamic  equations  are  in  mode  coupled  form  and  are  obtained 
by  expanding  the  surface  displacement  and  velocity  potential 
in  Fourier  series  (Hasselmann,  1961,  1963;  West,  et  al.,  1974). 
Our  treatment  begins,  as  does  that  of  Hasselmann  (1961,  1963) , 
with  the  time  dependent  equations  for  these  Fourier  coef¬ 
ficients.  Using  the  definition  of  a  local  power  spectrum 
due  to  Wigner  (1932) ,  a  transport  equation  is  obtained  that 
has  interaction  terms  of  lower  order  than  those  obtained 
by  Hasselmann  (1961,  1963). 


II .  The  Eigenmode  Equations  for  Surface  Waves 

The  difficulty  with  a  strictly  phenomenological  construc¬ 
tion  of  a  spectral  transport  equation  (Thomson  and  West,  1973) 
is  the  inherent  uncertainty  in  modeling  the  individual  inter¬ 
active  mechanisms.  Also,  because  of  the  nonlinear  nature  of  the 
interaction  process,  superposition  of  separate  mechanisms 
is  not  formally  justified,  although  it  may  be  pragmatically. 

To  obviate  these  complications  we  ure  the  dynamic  equations 
for  the  interaction  of  surface  waves  to  construct  a 
transport  equation.  We  are  concerned  with  a  description 
of  the  evolving  spectrum  of  surface  gravity  waves  due  both 
to  the  nonlinear  interactions  among  these  waves  and  to  the 
interaction  with  a  prescribed  oceanic  current.  Such  a 
current  might  represent  tidal  currents  at  the  mouth  of 
a  bay  or  estuary,  perhaps  wind -driven  currents,  etc.  We 
assume  the  characteristic  distances  over  which  this  current 
varies,  as  well  as  the  ocean  depth,  to  be  very  large  com¬ 
pared  to  the  wavelength  of  the  surface  waves  being  studied. 

The  basis  for  our  development  will  be  a  set  of 
eigenmode  equations  similar  to  ones  published  previously  (West, 
et  al. ,  1974).  Incompressible,  irrotational  flow  is  assumed. 

The  fluid  velocity  is  thus  expressed  as  the  gradient  of  a 
velocity  potential  <J>, 

<J>  =  <t>  +  $  .  (1) 
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In  Eq.  (1)  the  velocity  potential  is  given  by  a  linear 
superposition  of  the  short  wavelength,  high  frequency 
surface  waves  represented  by  <j>  and  the  long  wavelength, 

/v 

slowly  varying  prescribed  current  represented  by 

The  undisturbed  ocean  surface  is  assumed  to  coincide 
with  the  plane  z  =  0  of  a  rectangular  coordinate  system. 

The  z  axis  is  directed  upward  and  the  (x, y)  plane  lies  in 
this  surface.  The  horizontal  flow  associated  with  the 
prescribed  current  is 

A 

Jj(x,t)  =  Vs$  ,  z  =  0  ,  (2) 

where  Vs  is  the  gradient  operator  acting  in  the  horizontal 
plane,  x.  -  (x,y)  is  a  vector  in  this  plane,  and  we  assume 

/v 

<J>  to  vary  slowly  with  z  so  we  can  evaluate  JJ,  at  z  =  0. 

The  vertical  flow  of  the  current  ||  is  assumed  to  be  very 
small  compared  with  |u|  ,  and  the  surface  displacement 

H(x,t)  due  to  this  flow  is  also  assumed  to  be  correspon¬ 
dingly  small. 

When  surface  waves  are  present  the  equation  of  the 
sea  surface  is  of  the  form, 

z  =  H (x , t)  +  s (x,t)  (3) 
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where  £(x,t)  represents  the  short  wavelength  vertical  dis¬ 
placement  due  to  surface  gravity  waves.  Because  H(x,t) 
is  considered  to  be  very  small,  in  our  discussion  we  shall 
replace  Eq.  (3)  by  the  equation 

z  =  C  (x,t)  .  (4) 

(The  self-consistent  hydrodynamic  implications  of  the  assump- 
tions  made  about  H(x,t)  and  ^  are  discussed  in  Appendix  A.) 

In  West,  et  al.,  (1974)  the  velocity  potential  $  and  vertical 
displacement  z,  were  represented  as  discrete  Fourier  series 
in  a  rectangular  ocean  of  large  area  A Time-dependent  equations 
were  obtained  for  these  Fourier  coefficients,  similar  to  those 
previously  obtained  by  other  authors  (Phillips,  1960;  Benney, 

1962;  Hasselmann,  1961,  1963),  and  numerical  integration  of 
these  equations  was  described.  For  our  present  application, 
it  is  desirable  to  introduce  a  somewhat  modified  modal  analy¬ 
sis.  The  velocity  potential 

E  <Mx» z,t)  ,  at  z  =  ;(x,t;  , 

is  used  instead  of  <t>Q  ,  defined  as  <J>  on  the  plane  z  =  0. 

The  difficulty  with  the  use  of  is  that  the  shorter  wave¬ 
lengths  may  undergo  many  e-foldings  of  attenuation  between 
the  true  ocean  surface  and  the  plane  z  =  0. 
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We  represent  the  flow  by  the  complex  amplitude  Z(x,t) , 
defined  by  the  equations 


A  =  V  (Z+Z  )/2  , 

D  A 

;  =  i(Z-Z*)/2  .  (5) 


Here  Vx  is  the  "velocity  operator"  (g  is  the  acceleration 
of  gravity) 


vx  =  (g/CD)*5  =  o)x/(D 

<d  =  .  (6) 

b 


These  quantities  are  assumed  to  operate  on  functions  ex¬ 
pressed  as  Fourier  series,  for  which  the  proper  operation 
is  self-evident.  For  example,  we  assume  Z  to  be  a  function 
defined  in  a  rectangular  ocean  of  area  Aq  and  write 


Zik  "jx 

A  (JO  e 

k 


(7) 


where  the  Fourier  coefficients  A  (Jc)  are  time  dependent. 
Thus , 


VxZ 


Zik  *x 
vk  e  ~~ 


t 
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where 


vk  =  (g/k)*5  =  wk/k 


is  the  phase  velocity  of  a  small  amplitude  surface  gravity 

wave  of  wavenumber  k  in  deep  water.  The  corresponding  angula 

ij 

frequency  is  cok  =  (gk)  . 

In  our  rectangular  two-dimensional  space  representing 
the  quiescent  ocean  surface,  the  Fourier  exponentials 
satisfy  the  relations 


d2x 


exp(ik*x) 
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e:<p(ik*x) 


S(x) 


(9) 


where  6  is  the  Kroneckerand  6  (x)  the  Dirac  delta  function 
k 

The  prescribed  current  is  given  the  Fourier  repre¬ 
sentation 


U (x,t) 


Y  U(K)  cosIK-k  -  nKt) 
K 


(10) 


where  Q  is  some  (presently)  unspecified  function  of|K|and 
K 

the  vector  mode  amplitudes  U (K)  are  also  unspecified. 


It  is  straightforward  to  obtain  from  the  fluid  dynamic 
equations  the  time-dependent  differential  equations  satisfied 
by  the  Fourier  amplitudes  A(JO  of  Eq .  (7).  This  procedure  is 

outlined  in  Appendix  B  for  the  case  that  there  is  no  prescribed 
current  or  U  =  0.  The  additional  terms  required  to  account 
for  this  current  are  obtained  in  Appendix  A  in  an  approxi¬ 
mation  that  keeps  only  terms  linear  in  the  surface  wave 
amplitudes.  The  resulting  equations  are  [here  A  =  dA/dt] 


A(k)  +  ia)k  A  (JO  =  TW(A)  +  Ty  (A)  +  T2(A)  +  T3  (A)  +  ....  (ID 

In  obtaining  Eqs .  (11)  we  have  neglected  surface  tension 

and  have  supposed  the  ocean  to  be  much  deeper  than  the  longest 

wavelengths  considered.  We  must  therefore  set  equal  to 

zero  those  amplitudes  corresponding  to  capillary  waves  or 

corresponding  to  wavelengths  comparable  to  or  Greater  than 

the  depth,  e.g.,  in  Eq.  (11)  ,  T  models  the  effect  of  wind 

w 

and  viscosity.  Based  on  a  model  of  Miles  (1957,  1950)  (see, 
also,  Phillips,  1969),  we  adopt  for  this  the  simple  linear  ex¬ 
pression  fan  explicit  derivation  was  given  in  West,  et  al.  (1974) 


TW(A)  -  [.a*JS/(2V  "  vk2]A(^) 


(12) 


Here  a  is  a  vector  parallel  to  the  wind  direction  and  having 
a  magnitude  dependent  on  the  wind  speed. 
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The  quantity  T„  describes  the  coupling  to  the  pre¬ 
scribed  current.  As  obtained  in  Appendix  A,  this  has  the 

form 


T 


U 


C(-)  (k,K,)  AQc-K)  exp(-i«Kt) 

K 

+  C(  +  )(k,K)  h&+V  exP  ( i^Kt } 


(13) 


where 


(t)(k,K)  =  Is  0(K)-[i  +  Os.±K)  u)|k;K-J 


(14) 


The  term  T2  in  Eq.  (11)  represents  the  nonlinear  wave-wave 
interaction  of  &  (A2 )  .  This  is  derived  in  Appendix  B  and 
has  the  form 


t2<A> 


All)  Alp)  +  rf'  ^Alt,)  A*<-£> 
(-H)  A*  (-p)  1  •  (15) 


The  explicit  expressions  for  the  coefficients  r  are  given 

in  Eq.  (B. 7)  of  Appendix  B. 

Finally,  the  term  ri’3  describes  nonlinear  wave  inter 

actions  of  #(A3) .  This  is  also  shown  in  Appendix  B  to 
have  the  form 


'  3 


Vw-jg  A(^  A(£*  A  .te> 


r7,P,-n  A<i)  *<£>  A<-JJ)  +  if  ~'~A<1>  A*  ( — p )  A*  (n) 


+  r  a  (-£)  a  ( 


-p)  A* (n) j 


(16) 


k  n 

Of  the  above  coe cficients,  only  r'^wiii  J  . 

y  \&,p  11  be  needed  m  this 

Paper.  This  is  given  i„  Eq.  (B.8)  of  Appendix  B.  The 
remaining  terms  in  Eg.  (16)  tend  to  have  rapidly 

exponentials  and  are  not  expected  to  contribute  signifi¬ 
cantly  to  transfer  of  excitation  between  modes  in  Eg.  ,11,  . 

The  final  dots  in  Eq.  (U>  indicate  that  we  have 
included  only  terms  to  ©(A3)  in  the  interaction.  That  is, 
we  suppose  the  amplitudes  to  be  sufficiently  small  that 
terms  of&|A4|  and  higher  may  be  neglected. 


III.  Correlation  Functions  and  the  Power  Spectrum  of  the 
Wave  Amplitude 

We  wish  to  construct  spectral  transport  equations 
from  Eq.  (11).  To  do  this  we  first  introduce  an  ensemble 
average,  indicated  as 

”<•••>"  - 


r 


f 


i 


r 


over  many  observations  of  the  sea  state.  In  this  way  we 
can  construct  a  hierarchy  of  average  quantities  such  as 

<AQc)>  *  <  A  (JO  A(£)>  ,  etc. 

Using  Eq.  (11) ,  we  then  obtain  a  corresponding  hierarchy  of 
coupled  equations  for  these  quantities.  To  close  this  set- 
we  need  a  statistical  postulate  to  express  the  higher  order 
correlation  functions  in  terms  of  lower  order  correlation 
functions. 

The  postulate  which  we  adopt  is  that 

<AU)  A  (jp)  A*(k)  A*(n)>  =<AU)  A*(k)><A(p)  A*(n)> 

+  <AU)  A*  (n)>  <  A  (p)  A*  (k)  >  ,  (17) 

^  ^  ^ 

with  all  other  fourth-order  correlation  functions  vanishing. 
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X 


We  shall  also  neglect  as  small  all  correlation  products 

* 

involving  more  than  four  factors  of  A  and  A  . 

The  above  oostulate  permits  us  to  use  Eq.  (11) 

to  express  the  averages  of  products  of  one,  two,  and  three 

*  . 

A's  in  terms  of  the  v  A  (_£)  A  (n.)  >  and  to  also  obtain  an 
equation  to  determine  this  second  order  correlation 
function. 

The  postulate  (17)  has  arisen  in  a  variety  of  guises 
in  other  applications  of  statistical  mechanics,  perhaps 
the  earliest  being  Boltzmann's  assumption  of  "molecular 
chaos".  It  is  clearly  only  an  approximation.  In  the 
remainder  of  1.  iis  paper  we  shall  accept  this  postulate  and 
shall  not  attempt  to  assess  its  validity  here. 

At  this  point  it  is  convenient  (but  not  necessary) 
to  remove  the  term  T2  from  Eq.  (11)  with  a  transformation 
on  the  Fourier  amplitudes.  We  write 

A(J<)  =  a(k)  +  G(k)  ,  (18) 

where  G  satisfies  the  equation 

G(k)  +  iu)k  G (JO  =  T2(a)  .  (19) 

This  equation  may  be  formally  integrated  to  give 
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(k)  =  e 


k  tj  /V 

Thus  G  is  of  ).  The  difference 


T2(a)  dt' 


(20) 


T^(a)  =  T2(A)  -  T2(a)  # 

expressed  as  a  functional  of  the  a's  is  of  #(a3) .  This 
lets  us  finally  re-write  Eq.  (11)  in  the  form 

a(Jk)  +  iwk  a(k)  =  Tw(a)  +  T(J(a)  +  T^(a) 

where 

T3  (a)  =  T3(a)  +  T'(a) 

and  we  have  dropped  terms  of  order  higher  than  #(a3) .  We 
have  also  dropped  the  higher  order  terms  in  Tw  and 
which  is  consistent  with  our  use  of  only  simple  linear 
models  for  these. 

Equation  (21)  contains  only  terms  with  an  odd 
number  of  factors  of  the  a's.  Thus,  it  is  consistent  with 
this  equation  to  require  that  the  average  of  any  product 

with  an  odd  number  of  a-f actors  vanish.  We  have  then,  for 
example, 

<a(k)>  =<a(k)  aU)  a*  (n)>  =  0  ,  (23) 


,  (21) 

(22) 
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etc . 


The  corresponding  averages  of  the  A's,  obtained  with 
the  use  of  Eq.  (18),  ar-  not  expected  to  vanish.  Since 
we  are  neglecting  higher  than  fourth  order  correlation 
functions,  we  may  simply  replace  the  A's  by  the  quantities 

a  in  Eq.  (17) . 

We  shall  formally  suppose  that  Tw  ,  Ty  ,  ant  T3 
are  of  the  same  order  of  smallness  in  Eq.  (21).  This 
permits  us  to  evaluate  Eq .  (20)  in  a  simple  approximation, 

writing 

-iw,t 

a  (k)  -  q(k)  e 


and  considering  the  time  variation  of  the  q's  to  be  very 
slow.  Then,  we  obtain 


°  <i> 3 1 L  Vi-p 

I'Z 


r7p  a(^  a(£} 

«*>» +  %  - 


(l)  a  (-p)  iS '  '  ~i?  * 

y  a  I  •• 


^  “  “p  -  wk 


“j,  +  up  +  “k 


a  (-£)  a  (-JRL 


(24) 


Substitution  into  Eq. 


(22)  lets  us  finally  write 


(a) 


k,n 

6k+n-Si-p  C^fP 


a  (A)  a  (p)  a*  In) 

^  ^ 

+  terms  not  needed. 


(25) 
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The  "terms  not  needed"  here  have  the  form  of  the  final  three 

terms  in  Eq.  (16).  They  have  rapidly  oscillating  exponentials 

and  will  not  contribute  to  the  calculation  of  the  correlation 

functions  <a(U  a*(k)>  ,  according  to  our  statistical  postu- 

late  (17)  and  its  accompanying  postulates.  Were  we  to 

evaluate  <a(£)  a (p)  >  ,  on  the  other  hand,  we  would  require 

k  *  n 

some  of  these  other  terms  in  T^.  The  coefficient  is 

giver  in  Appendix  C. 

Following  Wigner  (1932)  we  now  introduce  t*^  power 
spectrum  of  the  a's  with  the  definition 


F  (x  ,Jk)  =  ^  ^  e1^*-  <a  (k+^/2 )  a*  (k-p/2 )  > 

£ 

=  (2AQf1Jd2r  e~1~£  <j(x+r/2)  ^  (*-£/2)>  .  (26) 


Here 


1 


E a(- 


(x,t)  =  )  a(k)  e 


(27) 


2 

which  dif feres  from  the  quantity  (7)  by  terms  of  0  'a  ) 
Using  Eq.  (5),  we  see  that 


^  Ffx,J$)  =  <S2  (£)>  +  terms  of  #(a4)  . 

& 


(20 
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In  the  next  section  we  shall  use  Eq.  (21)  and  the 
statistical  postulate  (17)  to  obtain  an  equation  for 
F(x.k).  The  other  correlation  functions,  <a(^)  a  and 

its  complex  conjugate,  may  then  be  evaluated  in  terms  of 
the  F's  [of  ^(F2)]  by  an  argument  similar  to  that  by  which 
Eq.  (24)  was  obtained.  Using  the  above,  power  spectra  for 
wave  energy  or  wave  amplitude  can  be  constructed.  We  may 
then  consider  F(x,k)  to  represent  an  approximation  to  tne 
power  spectrum  for  wave  amplitude.  The  precise  power 
spectrum  for  wave  amplitude  will  contain  additional  terms 
of  0tF2),  whi-h  can  be  readily  evaluated  as  just  described. 

For  most  applications  it  is  convenient  to  change  from 
discrete  to  continuum  norma li at ion  by  replacing  the  sum 
over  discrete  wavenumbers  by  integrals  with  the  substitution. 
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in  the  approximation  of  Eq.  (28).  The  function  4>(x,k) 
defined  by  Eq.  (30)  provides  a  generalization  of  the  spectral 
function  Y(k)  described  by  Phillips  (1969)  for  a  spatially 
homogeneous  ocean. 

In  practice,  the  Wigner  spectral  function  [Eq.  (3o)j 
is  useful  only  if  f  (x,Jc)  varies  very  slowly  over  distances 
comparable  to  k  for  all  k  of  interest.  For  oceanic 
applications  this  condition  is  usually  well  satisfied 
except  near  physical  discontinuities  (such  as  the  shore). 

We  thus  introduce  a  characteristic  distance  W  over  which 
H'Cgjjc)  varies  appreciably  and  assume  that 

k  >>  W_1  (32) 

for  those k  of  interest  .  Referring  back  to  Eq.  (26)  we 
see  that  Eq.  (32)  implies  that 

<^a(k+£/2)  a*  (k-;p/2^>«j  0  (33) 

for  |jd  |  >>  W  x. 

The  spectrum  of  energy  per  unit  area  is,  correct  to 
second  order  in  the  a's, 


E  (x,k)  =  p  g  vF(x,k) 

where  pQ  is  the  sea  water  density. 

—  - 

This,  for  example,  implies  that  k  >>  K  in  Eq.  (13). 
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If  we  write  aUs.t)  to  indicate  the  explicit  time  depen. 
dence  of  the  a(k)'s,  we  may  express  the  spectral  distribution 
for  wavenumber  and  frequency  in  the  form 


T  (x.k.t.w)  =  |ft0/[2(21I)3jJ^Jd 


dt  exp  (ip* x,  +  ) 


x  ^a(k+p/2,  t+x/2)  a*(k-£/2,  t-x/2>^> 


The  equation  satisfied  by  this  quantity  is  more  complicated 
than  that  derived  in  the  next  section  for  4-(x,k)  and  will 
not  be  described  here. 


IV.  The  Spectral  Transport  Equation 


In  the  preceeding  section  we  introduced  the 
spectral  function  in  terms  of  an  ensemble  average 
over  products  of  the  eigenmode  amplitudes.  In  Appendix  A 
we  presented  the  dynamic  equations  for  the  mutual  inter¬ 
action  of  these  surface  eigenmodes  and  their  interaction 
with  a  prescribed  surface  current.  In  this  section  we 
synthesize  these  approaches  to  construct  an  equation  for  the 
evolution  of  the  spectral  function  To  obtain 

this  equation  we  differentiate  the  first  form  of  Eq.  (26) 
with  respect  to  time; 

(S,jc) 

3t 


The  time  derivative  of  the  complex  amplitude  a  (JO  can  be 
eliminated  from  Eq.  (34)  by  substitution  from  Eq.  (21). 
We  then  obtain  on  the  right-hand  sidt  of  Eq.  (34)  a  sum 
of  terms  involving  correlation  functions  such  as 

<aQs>  (35a) 


=  *5 


^  <J£+.£/2  }  a  *  (k-£/2  £>  +  <a  (k+p/2 )  a*(K  -£/  2 


exp (i£>x) 


(34) 
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and 


(jd)  a(l)  a  (n)  a  (k)^ 


(35b) 


The  spectral  function  F(x,JO  is  of  course  just  a  Fourier 
transformed  version  of  Eq.  (35a)  . 

The  classification  of  correlation  functions  given 
by  Eq.  (35)  suggest  that  we  rewrite  Eq.  (34)  in  the  form 


3  F  (j&  i  )  _  m  .  m 
3t  a  b 


(36) 


We  consider  first  the  rather  simple  T  term.  This  is 

a 


Ta  " 


-ill- 


-  U) 


|k-£/2|  +  “(iS+£/2>  +  “  <iS-£/2>] 


x  ^a(k+£/2)  a  (£-£/2)^ 


+  )  |c  ("([Jk+p/2,K)  (a  (k+£/2-K)  a*yc-£/2)^>  exp(-i^Kt) 

K  L 

+  C(  +  U+p/2,j<)  (a  (k+£/2+K)  a*  exp(iQRt) 

-  C  ("\k-£/2,K)  \a(£+p/2)  a*  Qc-^p/2-^)^  exp(if2Rt) 


-  C  (+^jk-p/2f<K)  <^a(k+p/2)  a*  (J^-^p /2+K exp(-ifiRt) 


exp  (ijj* 


(37) 


19 


where  a  (jk)  =i  |j**j£/2vk  “  vk  ]  • 

Because  of  the  variation  in  'V  (x,k)  indicated  by 

Eq.  (33)  and  the  assumption  that  |jc|  >>  |j<|  t  we  need 
keep  only  the  lowest  order  non-vanishing  terms  on  expanding 
the  functions  in  Eq.  (37)  in  £  and  K.  With  some  little 
algebra,  then,  we  find  that 


T  =  -  C,  'V  F(:,k)  +  (ofjc/v.  -  2vk  )  F(x,k) 
a  —*k  i  1 ~  - K 


-u.yx  f 


(x,k)  +  [vx(k-U)l  -|\  F(x,k)j 

-  ^Vx  *[us^’ Ji/  (  2  k  2  )]^  F  (*  rJO 


(38) 


where  £k  s  VkoR  ,  Vx  is  the  gradient  operator  in  the  hori¬ 
zontal  plane  and  Vk  is  the  corresponding  wave  vector  gradient, 

The  assumed  form  [  Eq .  (17)]  of  the  fourth-order 

correlation  function  allows  us  to  express  in  Eq.  (j6)  in 
terms  of  the  spectral  density  function  F  (k,jjc)  .  A  straight¬ 
forward  evaluation  leads  to  the  expression 


=  ~  ^  jd2y  F(y,L)  f i2r  ^  exp|i£-  (x-r)  +  iq*  (r-y)J 
t.  J  J  p ,  q 


k+p/2,L-g>/2  *'J±/2’k+%/2  ... 

*  ^/2,^/2-aFC£'^/2)  FC£'*+2/2) 


(39) 


-1 


We  expect  p and  q  to  have  magnitudes  of  order  W  ,  and  thus 
to  be  very  small  compared  with  k  and  L  in  Eq.  (39) .  We 
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can  therefore  expand  G>  and  F  in  terms  of  and  <3,  keeping 
only  first-order  terms.  We  detine^5^  and  2  by  the  equations 


p.  £,L 

~1  --  ’k«L,Ji 


q  =  1  im 
q+0 


k,L+3/2  k,L-q/2 

^L-q/2,k^  "  ^L+q/2'^"^ 


(40) 


and  find  that 


T. 


V"  f  JS'L  r  1 

1  r  1 

mL  Kfc* 

•  hk  F  ‘i£  - i  > 

L  l  L 

1  L  ^  J 

+  fi-,x[F(iS'i!  F&-£)]j 

Using  the  forms  of  Tfl 

and  given  by 

(38)  and  (41)  in  Eq .  (36) 

,  we  can  express  the 

(41) 


time  derivative  of  F(x,k)  in  terms  of  F  itself.  K  more 
convenient  expression  may  be  constructed  by  using  T  (jc,Jt)  however . 
Employing  Eqs.  (29)  and  (30)  and  with  a  little  re-arranging 
we  obtain, 


L  +  .  v  +  .  y 

it  dt  x  at  k 

In  Eq.  (42) 


YO^k)  =  S(k)  H*  (x ,  Jc )  +  (a*k/vk-2vkz)4' (x  ’J  . 

(42) 


s*  = 

k  • 


dt 


(43) 


&  =  -  v  H 

dt  x 
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with  the  boundary  conditions 


x  =  x0/  iS  =  ii0  at  t  =  0 

time  t  =  0 ,  let  us  assume  that 

4*  =  ^(Xfk) 

Then,  at  time  t 

»US,k>  *  '*0[20<&-l'U  •  JSoU-i.t)]  '  <47) 

where  Eq.  (46)  has  been  inverted  to  express  jcQ  and  JcQ 

as  functions  of  jc,  Js»  and  t. 

The  second  term  on  the  right  in  Eq.  (42)  represents, 

in  a  fairly  obvious  way,  the  implications  of  our  modelling 
of  wind  and  damping  forces.  The  first  term  in  the  function 
S  (Jc)  can  be  rewritten  in  tensor  notation  as 

^  r  1  f9U •  9Ui  \ 

-  2^  kikj/(2k  lalTT  +  3x7  j 

i ,  j  =  l  L  ' 

This  is  seen  to  correspond  to  the  radiation  stress  term  intro¬ 
duced  by  Longuet-Higgins  and  Stewart  (1960,  1961).  The  re¬ 
maining  portion  of  YS  represents  a  kind  of  "stress"  associated 
with  the  interaction  of  waves  of  wavenumber  k  with  the 
entire  spectrum. 
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Had  we  kept  higher  order  terms  in  Eq.  U*> 
ponding  hi,her  order  ter„s  would  have  been  -rained  in 
For  example,  fifth-order  terms  in 
have  led  us  to  third-order  terms  .similar  to  those  Found 
bv  Hasseltann,  19.1.  19«3>  i»  «•  ««).  To  what  order 


one  can  continue 
order  correlation 


and  yet  neglect  fifth  and  higher 
functions  is  presently  not  evident 


24 


v .  The  Phase  Velocity  of  a  "Test  Wave" 

The  group  velocity  at  wavenumber  Jc  is  obtained  from 

Eqs .  (43)  and  (44)  as 


& 

dt 


U  +  £i 


/- 

-  d2L  S1  4<C£,L) 


(48) 


The  first  two  terms  are  obvious.  The  third  term  represents 
the  influence  of  the  nonlinear  wave  interactions.  Since 
V(x,L>  will  in  general  be  asymmetric  due  to  the  influence 
of  wind  and/or  obstructions,  the  group  velocity  can  have 

a  component  not  parallel  to  k. 

Some  insight  can  be  obtained  into  Eq.  (48)  by 
considering  a  "test  wave"  interacting  with  a  spectrum  of 
ocean  waves  in  a  uniform  ocean.  We  imagine  the  test  wave 
to  be  mechanically  generated  with  identical  characteristics 
for  each  of  a  sequence  of  observations.  Thus,  we  write 


aU)  “  ao(ii)  +  a'  6k -l  ' 


(49) 


where  a0  is  a  random  variable  describing  the  ambient  sea 
and  a'  represents  the  small  amplitude  "test  wave".  We 
substitute  (49)  into  Eq.  (21),  neglect  the  prescribed 
current,  wind  and  viscosity  terms,  and  obtain  a  linear 
equation  for  a' <k) .  Because  we  have  assumed  a  uniform  ocean, 
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<a0li>  ao(-E’>  *  st-p  <!  a0  (i> 1 2  > 


Integration  with  respect  to  time  then  gives  us  the  angular 
frequency 


(50) 


where  we  have  indicated  no  ^-dependence  of  41 .  The  phase 
velocity  is  wk/k  and  the  group  velocity  deduced  from 
F,q.  (50)  is  in  agreement  with  that  of  Eq.  (48). 

To  illustrate  the  implications  of  Eq.  (50)  we  consider 
the  spectrum  of  Tyler,  et  al.  (1974),  which  is  based 
on  a  representation  proposed  by  Longuet-Higgins ,  Eartwright, 
and  Smith  (1963).  This  is 


Y(L)  S  (0.4xl0"2/L4)  [g(B)/n]  ,  for  kQ  <  L  <  k p  , 


=  0  for  L  <  kQ  or  L  >  kp  .  (51) 


Here  the  angular  variation  of  the  spectrum  is  given  by 
G  (B)  =  a  +  (1-a)  coss(L)(B/2' 


and 


dB 
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_  2 

In  these  equations  a  is  very  small  (^-10  ),  and  kQ  and  k 

are  the  respective  long  and  short  wavelength  cut-offs  of 
the  spectrum,  and  B  is  the  angle  between  L  and  the  wind 
direction.  Finally,  s (L)  is  a  function  of  wavenumber  which 
is  near  unity  for  short  wavelengths  and  becomes  quite 
large  compared  with  unity  near  L  =  kQ  . 

We  shall  evaluate  Eq.  (50)  for  wavelengths  shorter 
than  the  cut-off,  or 

k  >>  kQ  (52) 

In  this  case  the  principal  contribution  to  the  integral 

in  Eq.  (50)  comes  from  L-values  near  L  =  k  and  a  simple 

analytic  evaluation  is  possible. 

k,L 

The  coefficient  ~  is  obtained  from  Eq.  (c.l)  of 
Appendix  C.  For  k  >>  L,  this  is 

tfL,k  S  '  (t)L  VLJS-i  ’  <53) 

On  evaluating  the  integral  we  find  that 

uk  S  wk  [l  +  1.4  xlO-2  cosB’  (k/k^J  ,  (54) 

where  3  is  the  angle  between  k  and  the  wind  direction.  The 
group  velocity  obtained  from  Eq.  (54)  is 
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VI .  Wave  Shadowing  by  an  Island 

During  the  series  of  experiments  reported  in  Tyler,  et  al.^ 
(1974),  the  "shadow"  of  an  island  for  receding  waves  was  observed 
At  sufficiently  large  distances  from  the  island  this  shadow 
is  absent.  There  are  evidently  several  possible  causes  for 
the  filling  in  of  the  spectrum  away  from  the  island.  One 
of  these  is  nonlinear  wave  interactions,  which  we  now 

discuss  as  an  application  of  Eq.  (42). 

We  calculate  the  filling  in  of  the  spectrum  of 

wavenumbers  k  directed  aw,  y  from  the  island  and  in  its 
shadow.  That  is,  ¥  (x,k)  will  be  very  small  where  effec¬ 
tive  shadowing  occurs.  On  the  other  hand,  we  assume 
that  V (x, L)  =  ¥ (L)  will  not  have  much  x-dependence  for 
those  waves  JL  which  have  "missed"  the  island.  If  the 
shadowing  angle  is  small,  we  can  take  (we  now  suppose  that 
^  =  0  and  the  effects  of  wind  and  viscosity  can  be 

neglected) 


at  o  ~  0 

—  U '  At — 


in  Eq.  (42)  .  Equation  (48)  gives  the  group  velocity  with  which 
waves  of  wavenumber  k  propagate  into  the  shadow 


C(k>  -  Ck  -  |d2I,g1  YIL)  . 


*We  are  indebted  to  Professor  Walter  Munk  for  describing 
these  observations  to  us  prior  to  publication. 
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If  this  were  a  time-dependent  problem,  with  a  sharply 
outlined  shadow  at  t  =  0,  say  V  =  VQ(x,k),  then  at  time  t 
Eq.  (47'  would  imply  that 

Y(x,k)  =  fQ(x  -  C(k)t,  k)  .  (57) 

The  expression  (55)  would  lead  us  to  expect  a 
triangular  shadow  of  half  angle 

0i  =  1 . 4X10~2  sin  (28)  (k/k^1*  .  (58) 

When  waves  travelling  parallel  to  the  wind  are  shadowed 
by  the  island,  then  the  filling  in  of  the  spectrum  wi.1 1 
be  modified.  Should  a  significant  portion  of  the  spectrum 
be  in  the  island  shadow,  then  'ML)  in  Eq.  (56)  must  be 
appropriately  modified. 
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VII.  Other  Correlation  Functions 

The  correlation  function  < a  (£)  a  (k)  >  is  easily  ob¬ 
tained  in  terms  of  F,  or  equivalently,  < a (l)  a  (p) >  • 
Using  Eq.  (21)  (and  ignoring  TyJ  and  V  ,  we  obtain 

1  V  [  ~~~ 

I  *  +  i(u)k+to  )  <a(q)  a(k)>  =  ^ 

L  8.,p.n 

**+  *** 


+  6q+n+p-S-  CJL 

<4* 


x  <a(£)  a  CL)  a*  (jg)  a*  (n)  > 

q,p,n  *  *  .  , v  i 

~  ~  <  a  ( k )  a  (£)  a  (p )  a  CC.)  >  J 


=  T(k,5) 


Thus 


<a(q)  a(k)  >S  -i(wR+wq)  T(k,q)  • 


Use  of  Eq 
quantity . 


(17)  then  permits  explicit  evaluation  of  this 
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APPENDIX  A:  The  Interaction  of  Surface  Waves  with  a  Current 

In  this  Appendix  we  obtain  the  terms  in  Eq .  (11) 

which  represent  the  interaction  of  surface  waves  with  the 
prescribed  current  U.  Bernoulli's  equation  and  the  kine¬ 
matic  boundary  condition  at  the  surface  are,  respectively, 

+  H  (V$)2  +  gh  =  0  , 

a  t 


where  both  equations  are  evaluated  on  the  surface 


z  =  h  (x ,  t )  =  H  +  C 


(A. 2) 


If  we  extract  the  long  wavelength,  low  frequency 
part  of  Eq.  (A.l)  in  linearized  form,  there  results 


3  H  d± 

at  at 


+  gn  =  0 


(A. 3) 


both  evaluated  on  the  surface  z  0. 


A 


The  current  associated  with  $  is 

✓v  A  A  A 

U  =  V4>  =  U  k  +  U  / 


where  U  is  given  by  Eq.  (2).  We  have  assumed  that 


/A 


(A. 4) 


Our  conditions  also  imply  that 


|V  in  H|  <<  |Vs  &n  c| 


(A. 5) 


or  that  the  surface  curvature  due  to  the  current  is  much 

less  than  that  due  to  the  surface  waves. 

To  extract  the  high  frequency  linearized  terms  from 

(A . 1 )  we  note  that 


A. 


A 


2^ 

since  V  $  =  0. 
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Thus,  the  high  frequency  part  of  Eq.  (A.l),  linear 
in  the  surface  wave  amplitudes,  and  evaluated  on  the 
surface  z  =  H,  is 


at 


A 

+  u*  V4> 


+s 


3U 

IT  +  « 


0  , 


H 

at 


KT-  +  u 


V 


(V  <C)  *  (V  H)  +  5  (V  -U) 

S  b  b 


a^_ 
a  z 


(A. 6) 


Equation  (A. 6)  is  to  be  evaluated  by  replacing 
<J>(x,z,t)  by 


$  (x,H,  t)  =  <p  (x ,  t ) 

~  ti  ~ 


(A. 7) 


after  the  indicated  differentiations  are  performed. 

Following  a  treatment  of  Milder  (1973)  ,  we  re-express 

Eq.  (A.  6)  in  terms  of  <j>  .  For  example, 

H 


a 

at 


+ 


'z=H 


+  u*v 

—  s 


(A. 8) 


This  lets  us  replace  the  first  of  Eqs.  (A.  6)  by 


a 

at 


u*  v 

~  s 


<p  +  g 
H  y 


eff 


C  =  0 


(A. 9) 


au 

geff  =  g  +  TT 


(A. 10) 


35 


' 


/V 

The  unit  normal  n,  to  the  surface  is,  to  first  order  in  H 


n  =  k  -  V  H 

-w  -w  5 


(A. 11) 


Thus,  the  second  of  Eqs.  (A. 6)  is  equivalent  to  the 
equation 


Jt  +  £'Vs)  ?  +  ^Vs‘£)  =  «.V4» 


on  z  =  H 


(A. 12) 


To  express  this  in  terms  of  <)>  ,  we  use  the  relation 

2 

V  <J>  =  0  to  write 


n*  VcJ)  I 

1  z=H 


=  d>^H 


2 

3  H 


3x 


II 


and 

<D  = 


(A. 13) 


(A. 14) 


The  conditions  (A. 4)  and  (A. 5)  finally  give  us  the  equations 


(A  +-H‘Vs)  %  +  «  *  0  ' 

+  H.-7S]  5  +  C  7S-U  =®  *H  •  (A. 15) 
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APPENDIX  B:  The  Wave-Wave  Interaction 


In  this  Appendix  we  show  how  to  obtain  the  terms 
T2  and  in  Eq.  (11). 

Extracting  from  Eqs.  (A.l)  the  part  which  pertains 
to  surface  gravity  waves,  we  obtain  the  equations 


II  +  >5  (V<M2  +  gr.  =  0  ,  2  =  c 

§£ +  <v>-<v>  -  H  -  z  -  5  •  <Ba’ 

We  first  re-express  these  equations  in  terms  of  <J>(x,z,t) 

evaluated  on  the  surface  z  =  £(x,t);  that  is, 

4>s(x,t)  =  <}>[x,  5  (x,t)  ,  t]  .  (B.2) 


Then,  we  define 

w&'fc)  -  H|2=c 

and  re-write  Eqs.  (B.l)  as 

TT  "  W  H  +  ’s  (vs<t>s-wV)  +  h  w/  +  g;  =  0 

f£+ <Vs-MV>-V*  w  ■ 
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It  remains  to  express  W  in  terms  of  <j>  ,  which  is  a  special 

and  rather  simple  application  of  potential  theory  with  a  Dirichlet 
boundary  condition (Jackson ,  1962).  This  is  easily  done 
by  first  expressing  both  <f>s  and  W  as  Taylor  series  in  q 
about  the  plane  2=0.  Then  W  can  be  expressed  in  terms  of 
<j>  by  successive  substitution.  The  result  is 

W  =((D  4>s)-  [®(c®  4>s)“  (C®24>s)]  +  ^)[^®(^s)]  -  C^)2(C®<t>s) 

-  *s^lc2(®2<>s)]  "  ^((D34>s)j  •  <B*5> 

The  term  can  bj  eliminated  from  the  first  of 
3 1 

Eqs.  (B.4)  using  the  second  of  these  equations  and  W  eliminated 
from  both  using  Eq.  (B.5).  Finally,  a  first-order  equation 
for 

z  =  -U  +  V^1  <f>s  (B.6) 


can  be  obtained  by  differentiation  with  respect  to  time  and 
substituting  from  Eqs.  (B.4).  The  Fourier  expansion  [Eq.  (7) 
then  gives  us  the  terms  T2  and  T^  of  Eq.  (11). 

The  coefficients  in  T2  [Eq.  (15)J  are 
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ri7p  =  ff  [(WVk)Up+-i;  '  VkP-*\8>  -  vn 

Js»£  1  f  .1 

=  4  I ‘v£vp//vk)  "  vp(kp+~‘p)  +  vi  (kfc_iS‘i.)  * 

■&'*&'£  if  1 

r  =  8  I  (v£vp/vk}  Up+i/i>)  +  vp(kP+k-j>)  +  (ka+jc-^)  .  (b. 7 ) 

The  first  coefficient  in  Eq.  (16)  (the  only  one  required  in 
this  paper)  is 

n,  ,  r 

r£fp  =  T  [(u,n"u)p)  (k"li$riJ>  +  (a)n"U)£)  (k~|k-£l  ) 

"  (W^+Wp)  l-L+£l  (k  —  I  )  +  WpP(k~p)  +  w^Mk-Jl)  -  wnn(k-n) 

-  p-n_  -  a)  jt-n,  -  2wn  £*p  +  (a)jla)n/wk)k(n-|n-p|+J?.-|£+p|  ) 

+  (a)pU)n/a)k)k  (n-  ln-i.l+P-I^P+^1  )  -  (w£wp/u)k)k(p-|p-n|+«.-!<£->n| 

(B.  8) 

Since  the  condition  (32)  has  been  used  in  our  derivation  of 
Eq.  (42),  we  mast  restrict  ourselves  to  wavelengths  small 
compared  to  the  length  parameter  W.  To  do  this,  we  suppose 
the  coefficients  (B.7)  and  (B.8)  vanish  if  any  of  their  wave- 
number  arguments  violate  the  condition  (32) . 


APPENDIX  C:  The  Coefficients  in  Equation  (25) 

For  reference  we  quote  the  form  of  the  coefficients 
of  the  a’s  in  Eq.  (25) : 
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r*Js 

p-n,n 

k  2-n,n 

II 
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2  ,p-n 

pTW  ~ 
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~  ~ 

p,2-n  2 

A*  rSmf 

c 
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1 

a 

3 

_ l 

Q.I 
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CO  .  —  CO  —CO  I  „  I 

2  n  2-n 

*■%*  1 
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2 

-  t  '■'* 

n-p»P 

r  — 
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k,n-2 

r 

p 

n-2  2 

r~  ~  ~ 

n 

(D  -0J  -0J 

n  p 

|p-n| 

^  Ol)  -co  - 

n  2 

CO  1  1 

Urill 

k,n  ^2+p 

ri+jc  r_£  ,p  "  (i+p) 

0)  I  .  -t.  +  2  Co'-'+uT+CO  •  (C,1) 

|  p+2  I  A,  p  2,  p  Jl- +p 


For  the  evaluation  of  the  coefficients  (40)  certain  of  the 
terms  in  (C.l)  appear  to  be  singular,  corresponding  to  the 
resonant  excitation  of  arbitrarily  long  wavelengths.  In 
accordance  wich  the  discussion  following  Eq.  (B.8),  these 
terms  are  to  be  dropped ,  corresponding  to  the  assumed 
vanishing  of  the  T-coef f icients . 

For  the  evaluation  of  Eq.  (53) ,  one  should  note  the 
sequence  of  cancellations  of  the  terms  with  powers  of  k 
greater  than  the  first. 
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management  assistance  in  support  of  studies,  analyses, 
development  planning  activities ,  acquisition,  test, 
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